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Abstract. We obtain a class of parametric oscillation modes that we call K-modes with 
damping and absorption that are connected to the classical harmonic oscillator modes through 
the "supersymmetric" one-dimensional matrix procedure similar to relationships of the same 
Cf^ ' type between Dirac and Schrodinger equations in particle physics. When a single coupling 

parameter, denoted by K, is used, it characterizes both the damping and the dissipative features 
of these modes. Generalizations to several K parameters are also possible and lead to analytical 
■ results. If the problem is passed to the physical optics (and/or acoustics) context by switching 

I from the oscillator equation to the corresponding Helmholtz equation, one may hope to detect 

^2 . the K-modes as waveguide modes of specially designed waveguides and/or cavities. 

o 

g : PACS number(s): 12.60Jv, 1 l.SO.Pb arXiv number: math-ph/0402065] 

O '. 1. Introduction. 

"7^ I Factorizations of differential operators describing simple mechanical motion have been only 
PLi- occasionally used in the past, although in quantum mechanics the procedure led to a vast 
fH ! literature under the name of supersymmetric quantum mechanics initiated by a paper of Witten 
|[1 1. However, as shown by Rosu and Reyes |2|, for the damped Newtonian free oscillator the 
factorization method could generate interesting results even in an area settled more than three 
centuries ago. In the following, we apply some of the supersymmetric schemes to the basic 
classical harmonic oscillator. In particular, we show how a known connection in particle 
^ ■ physics between Dirac and Schrodinger equations could lead in the case of harmonic motion 
to chirped (i.e., time-dependent) frequency oscillator equations whose solutions are a class 
of oscillatory modes depending on one more parameter, denoted by K in this work, besides 
the natural circular frequency ujq. The parameter K characterizes both the damping and the 
losses of these "supersymmetric" partner modes. Moreover, we do not limit this study to one 
K parameter extending it to several such parameters still getting analytic results. Guided by 
mathematical equivalence, possible applications in several areas of physics are identified. 

2. Classical harmonic oscillator: The Riccati approach. 
The harmonic oscillator can be described by one of the simplest Riccati equation 

u' + + kujI = 0, K = ±l , (1) 

where the plus sign is for the normal case whereas the minus sign is for the up side down case. 
Indeed, employing u = one gets the harmonic oscillator differential equation 

w + kuJqW = , (2) 
with the solutions 

W+ cos(u;ot + (f^) if K = 1 



Wb 



W_sinh(t<;ot + (f-) if k = — 1 
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where W± and Lp± are amplitude and phase parameters, respectively, which can be ignored in 
the following. 

The particular Riccati solution of Eq. dl]) are 



u 



j — co'otan(co'ot) ifK = l 
P I u;ocoth(co'ot) if K = —1 



It is well known that the particular Riccati solutions enter as nonoperatorial part in the 
common factorizations of the second-order linear differential equations that are directly 
related to the Darboux isospectral transformations LIJ. 
Thus, for Eq. Q one gets (Dt = ^) 

(Dt + Up) (Dt - Up) w = w" + (-Up - uj)w = . (3) 

To fix the ideas, we shall use the terminology of Witten's supersymmetric quantum mechanics 
and call Eq. Q the bosonic equation. We stress here that the supersymmetric terminology 
is used in this paper only for convenience and should not be taken literally. Thus, the 
supersymmetric partner (or fermionic) equation of Eq. @ is obtained by reversing the 
factorization brackets 

(Dt - Up) (Dt + Up) Wf = w" + (Up - uj)w = w" + ^2(t)w = , (4) 
which is related to the fermionic Riccati equation 

u' - u^ - uj'j{i) = , (5) 
where the free term ujf is the following function of time 

2, ^ I 2 J ^o(^l ~ 2tan^u;ot) if k = 1 

(t) = Up - Up = j ^2(1 _ 2coth2cUot) if = -1 . 

The solutions (fermionic zero modes) of Eq. (|4|) are given by 



-^0 



if K 



Wf 



COS (a;ot) 



and thus present strong periodic singularities in the first case and just one singularity at the 
origin in the second case. These 'partner' oscillators, as well as those to be discussed in 
the following, are parametric oscillators, i.e., of time-dependent frequency. Moreover, their 
frequencies can become infinite (periodically). In general, signals of this type are known 
as chirps. 'Infinite' chirps could be produced, in principle, in very special astrophysical 
circumstances, e.g., close to black hole horizons yj. 

3. Matrix formulation. 

Using the Pauli matrices (^y = { ^- r, ) and = { i ) , we write the matrix equation 



V^W = [ayDt + a.(iup)]iy = , (6) 

where W = ^ ^ is a two component spinor. Eq. Q is equivalent to the following 

decoupled equations 

(iDt + iup)wi = (7) 

(-iDt + iup)w2 = . (8) 
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Solving these equations one gets wi oc cuq/ cos(a;ot) and W2 oc uo^ cos(co'ot) for the /t = 1 case 
and wi oc co'o/sinh(c<Jot) and W2 oc co'osinh(co'ot) for the k = —1 case. Thus, we obtain 



W 





(9) 



This shows that the matrix equation is equivalent to the two second-order linear differential 
equations of bosonic and fermionic type, Eq. Q and Eq. Q, respectively, a result quite well 
known in particle physics. Indeed, a comparison with the true Dirac equation with a Lorentz 
scalar potential S(x) 

[-i(TyD, + a,(m + S(x))]iy = EM/, (10) 

shows that Eq. Q corresponds to a Dirac spinor of 'zero mass' and 'zero energy' in an 
imaginary scalar 'potential' iup(t). We remind that a detailed discussion of the Dirac equation 
in the supersymmetric approach has been provided by Cooper et al Q in 1988. They showed 
that the Dirac equation with a Lorentz scalar potential is associated with a susy pair of 
Schroedinger Hamiltonians. This result has been used later by many authors in the particle 
physics context L6J. 



4. Extension through parameter K. 

We now come to the main issue of this work. Consider the slightly more general Dirac-like 
equation 

V^W = [(TyDt + a,(iup + K)]W = KW , (11) 

where K is a (not necessarily positive) real constant. On the left hand side of the equation, 
K stands as an (imaginary) mass parameter of the Dirac spinor, whereas on the right hand 
side it corresponds to the energy parameter. Thus, we have an equation equivalent to a Dirac 
equation for a spinor of mass iK at the fixed energy = iK. This equation can be written as 
the following system of coupled equations 

iDtWi + (iup + K)wi = Kw2 (12) 

-iDtW2 + (iUp + K)w2 = Kwi . (13) 

The decoupling can be achieved by applying the operator in Eq. (13) to Eq. (12) . For 
the fermionic spinor component one gets 



w 



1 -^0 



2 N 2K 
2 tan uiQt) + 1 — tan c^ot 

2cotn tuot) + i — cothcuot 



Wi 



Wi 



for K 



for K 



whereas the bosonic component fulfills 



. ■2K 

1 — 1 — tanc^ot 

2K , 

1 — i — cotncuot 



Wo 



for K = 1 



for K 



■1 . 



(14) 
(15) 

(16) 
(17) 



The solutions of the bosonic equations are expressed in terms of the Gauss 
hypergeometric functions 2F1 



^ Jp-|)Jq-|) p 

a+Zi Z2 2-r 1 



_^^e-2ip.4(p-i)^-(p-|)^(q-|)^p^ 



p + q,p + q- l,2p;-^zi 



q-p,q-p 



1,2 



2p;-2-i 



(18) 



r + s, r + s + 1, 1 + 2r; -Z3 
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and 

W2 (t; q;_, = a-z\-z\ 2F1 
+/3_4^Z3--z^2Fi 

where the variables (i = 1, 4) are given in the following form: 

Zi = itan(cJot) — 1, Z2 = itan(cjot) + 1, Z3 = coth(ci;ot) + 1, Z4 = coth(a;ot) — 1, 
respectively. The parameters are the following: 



s - r + l,s - r, 1 - 2r; -Z3 



(19) 



1 / 2K \ 1 / 2K \ 1 2K 1 2K 



The fermionic zero modes can be obtained as the inverse of the bosonic ones. Thus 

1 1 



(20) 



w^(t;a+,/3+) ' W2 (t; q;_, /?_) ■ 

A comparison of Wj*" with the common 1 / cos t fermionic mode is displayed in Figures (3) 
and (4). 

In the small K regime, K ^ 6i;o, one gets 



,(p-|)„(q-5) 



w+(t;a+,/3+)R^a+zr^^Z2 '^Fi 2 , 1 , 2 - -; --zi(t) 



K 1 



and 



W2(t;a_,/3_) «ia_z^z|2Fi 1, 2, 2 + i— ; -Z3(t) 



K 1 



(21) 



+ 



/3_4'-Z3-'^zl2Fi 



K K K 1 , , 
1 - 1—, -1—, -1— ; -Z3(t) 

LOq LOq Z 



(22) 



Examining the bosonic equations, one can immediately see that the resonant frequencies 
acquired resistive time-dependent losses whose relative strength is given by the parameter K. 
The fermionic equations having time-dependent real parts of the frequency can be interpreted 
as parametric oscillators which are also affected by losses through the imaginary part. 

5. More K parameters. 

A more general case in this scheme is to consider the following matrix Dirac-like equation 



i ] ^* + i 1 ] i iup + K2 



Wi 
W2 



k; 
k; 



Wi 
W2 



The system of coupled first-order differential equations will be now 

- iDt -\- iUp K2 W2 = K'^wi 
iDt + iUp Ki wi = K2W2 



(23) 

(24) 
(25) 
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and the equivalent second-order differential equations 



-iAK 



Dtw,+ ±DtUp+i(Ki+K2)up+(KiK2 - w, = ,(26) 



where the subindex i = 1,2 and AK = Ki — K2. Under the gauge transformation 

1 r n \ ] 

■ ■ " Z,;(t)e^ 



W,; 



Z^exp ( -- 



iAK 



dr 



,iitAK 



one gets 



Dt% + Q,(t)Z, = 0, 
where the 'potentials' have the form 



g,(t) = ±DtUp + i(Ki + K2)up + (KiK2-K;K 



n 1 






-iAK 


~4 





(27) 
(28) 
'(29) 



(5i,2 are functions that differ from the nonoperatorial parts in Eqs. (I30l-l33t only by constant 
terms. Indeed, one can obtain easily the following equations: 
For the fermionic spinor component one gets 

^2 KiK2-K;k; , .K1 + K2 



for K = 1, and 



. (K2-K/ 



+ 1- 



1 - 2coth^a;ot + 



2 ^ , (K2-K1)' , K1K2-K1K2 , .K1 + K2 



+ 



^0 



for K = —1. 

The bosonic component fulfills 



D?Z+ + ^^ 



K2 -Ki 2 
1+ ^ , 2 + 



KiK2-k;k; .K1 + K2 



+ 1- 



tancuot ZJ = 0(30) 



cothtuot = 0(31) 



^0 



UJQ 



for K = 1, and 



D?Z2 



UJr 



1 - 



(K2-Ki)2 KiK2-K;k; .K1 + K2 



^0 



Wo 



tan Uot 



coth cjot 



0,(32) 



0,(33) 



for K = —1. When Ki = K2 = K one gets the particular case studied in full detail above. 
The more general 'bosonic' modes have the form: 



Z^(t; a+,/5+) = «+[tan(co'ot) — ij^'-'o [tan(cjot) + ij^'-'o 
2F1 



, + 1,1 + — ; - tan Wot - 1 



[tan(u;ot) - ii 

^^2 — ^1 ^2 — ^1 



2F1 



4a;o 



4a;n 



1,1- 



'0 [tan(wot) + i] 

(tan(wotj — ij 



2a;o 2 



(34) 



and 



^3 04 

Z2 (t;a_,/3_) = a_[coth(wot) - 1] [coth(wot) + 1]*^ 



2F1 



—J, + 1, , 1 + -- coth Wot 

4iWo 4u;o iwo 2 



+/5_(-l) 2"o42"o [coth(wot) - 1] [coth(wot) + l]*"o 



4wo 



4wo 



1,1 



^;-^(coth(wot) - 1) 



(35) 
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where 



1^1 = [AujI + (Ki + K^f + 4[(Ki + K2)cuo - K'^K^ 
^2 = (4cUo + (Ki + K2)' - 4[(Ki + K2)cuo + K;k; 
= (AuI - (Ki + K^f - 4[i(Ki + K2)^o - KK 



^4 = U^l - (Ki + K^f + 4[i(Ki + K2)^o + K;k: 



1/2 
1/2 
1/2 
1/2 



6. Applications. 

6.1 Waveguides. 

In view of the correspondence between mechanics and optics, one can also provide an 
interpretation in terms of the Helmholtz optics for light propagation in waveguides of special 
profiles. The supersymmetry of the Helmholtz equation has been studied by Wolf and 
collaborators |7|. To get the waveguide application, one should switch from the temporal 
independent variable to a spatial variable t ^ x along which we consider the inhomogeneity 
of the fiber whereas the propagation of beams is along another supplementary spatial 
coordinate z. Thus, we turn the equations (14-17) into Helmholtz waveguide equations of 
the type (we take c = 1) 

[dl + dl+uy{y.)]^{^,z)=Q, (36) 

where the modes z) can be written in the form wi 2(x)e~^'^^^ for a fixed wavenumber k^ 
in the propagating coordinate that is common to both wavefunctions and the index profiles 
correspond to two pairs of bosonic-fermionic waveguides and are given by 

2K 2K 

n^(x) ~ 1 — i- — tan(kox) , n£(x) ~ — (1 + 2 tan^(co'ox)) — i- — tan(kox) ,(37) 

ko kg 

and 

2K 2K 

nb(x) 1 - i— coth(kox) , Yi}(yi) ~ 1 - 2coth2(kox) + i— coth(kox) , (38) 

ko ko 



respectively. In our units ko = cuo. Eqs. (0711381 can be obtained from Riccati equations of 
the type {c ^ 1) 

a;o%%(x)/c2 = k2^R.-R2, (39) 

where R(x) are Riccati solutions directly related to the Riccati solutions discussed in the 
previous sections. 

According to Chumakov and Wolf [7| a second waveguide interpretation is possible 
describing two different Gaussian beams, bosonic and fermionic, whose small difference in 
frequency is given in terms of a small parameter e (wavelength/beam width), propagating in 
the same waveguide. In this interpretation, the index profile is the same for both beams. 
For illustration, let us take the normal oscillator Riccati solution in the space variable x, 
i.e., tankox that we approximate to first order linear Taylor term kox. Then, the two beam 
interpretation leads to the following Riccati equation (for details, see the paper of Chumakov 
and Wolf) 

ujiy{^) - ^o'n'(O) = Tko - k2x2(l ^ e) . (40) 

An almost exact, up to nonlinear corrections of order and higher, supersymmetric pairing 
of the z wavenumbers (propagating constants) occurs, except for the 'ground state' one. As 
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noted by Chumakov and Wolf, supersymmetry connects in this case light beams of different 
frequencies but having the same wavelength in the propagation direction z. This approach is 
valid only in the paraxial approximation. Therefore, one should know the small x behaviour 
of the K-modes in order to hope to detect them through stable interference patterns along the 
waveguide axis. 

6.2 Cavity physics. 

Another very interesting application of the K-modes in a radial variable could be Schumann's 
resonances, i.e., the resonant frequencies of the sperical cavity provided by the Earth's 
surface and the ionosphere plasma layer [8|. The Schumann problem can be approached as a 
spherical Helmholtz equation [V^ + = with Robin type (mixed) boundary condition 

f^ls = C{uj)(j)s, where C{ijj) is expressed in terms of the skin depth 5 = ^J2/Jp^mv) of 
the conducting wall, fic is its permeability and a is its conductivity. The eigenfrequencies 
fulfilling such boundary conditions can be written as follows 



UJ 



(41) 



where / is a complicated expression in terms of skin depths and surface and volume integrals 



of Helmholtz solutions with Neumann boundary conditions ^ 



0. It is worth noting 



the similarity between these improved values of Schumann's eigenfrequencies and the K- 
eigenfrequencies. Moreover, using the Q parameter of the cavity, one can write Eq. (HTT) in 
the form 



1 



1 



1 



Q Q 



(42) 



This form shows that the modification of the real part of cu leads to a downward shift of the 
resonant frequencies, while the contribution to the imaginary component changes the rate of 
decay of the modes. 

We point out that Jackson mentions in his textbook that the near equality of the real and 
imaginary parts of the change in u'^ is a consequence of the employed boundary condition, 
which is appropriate for relatively good conductors. Thus, by changing the form of C{uj) 
that could result from different surface impedances, the relative magnitude of the real and 
imaginary parts of the change in cu'^ can be made different. It is this latter case that corresponds 
better to the K-modes. 

6.3 Crystal models. 

There is also a strong mathematical similarity between the K-modes and the solutions of 
Scarf's crystal model [9| based on the singular potential V{x) = — Vocosec^(7rx/a), where 
a is an arbitrary lattice parameter. For this model the one-dimensional Schrodinger equation 
has the form 



cosec 



{nx/a) 



^ = 



(43) 



For < X < a/2, the general solution is 



^=[f{xp^\F, 



- + -{s + X),- + -{s-Xy,l + s;f{x) 



+ 



l-l{s-X),\-^{s + Xy,l-s;f{x) 



(44) 



where f{x) = sin(7rx/a) corresponds to the Zi{t) functions, and s and A corresponding to 
-p and -q, respectively, are related to the potential amplitude and energy spectral parameter. 
Thus, by turning the K-oscillator equations into corresponding Schroedinger equations, one 
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could introduce another analytical crystal model with possible applications in photonics 
crystals. 

6.4 Cosmology. 

Two of the authors applied the K-mode approach to barotropic FRW cosmologies ifTOll . 
K- Hubble cosmological parameters have been introduced and expressed as logarithmic 
derivatives of the K-modes with respect to the conformal time. For K ^ the ordinary 
solutions of the common FRW barotropic fluids have been obtained. 

It is also worth noticing the analogy of the nonzero K oscillator case with the 
phenomenon of diffraction of atomic waves in imaginary crystals of light (crossed laser 
beams) ifTH . In fact, the K parameter is a counterpart of the modulation parameter Q 
introduced by Berry and O'Dell in their study of imaginary optical gratings. Roughly 
speaking, the nonzero K modes could ocurr in an imaginary crystal of time that could occur 
in some exotic astrophysical conditions. 

7. Conclusion. 

By a procedure involving the factorization connection between the Dirac-like equations and 
the simple second-order linear differential equations of harmonic oscillator type, a class 
of classical modes with a Dirac-like parameters describing their damping and absorption 
(dissipation) has been introduced in this work. While for zero values of the Dirac parameters 
the highly singular fermionic modes are decoupled from their normal bosonic harmonic 
modes, at nonzero values a coupling between the two types of modes is introduced at the level 
of the matrix equation. These interesting modes are given by the solutions of the Eqs. (l30b 
- and in a more general way by Eqs. (ITTl 34-35) of this work and are expressed in 
terms of hypergeometric functions. Several possible applications in different fields of physics 
are mentioned as well. Finally, similar to the fact that the PT quantum mechanics can be 
considered as a complex extension of standard quantum mechanics, we notice that what we 
have done here is a particular type of complex extension of the classical harmonic oscillator. 
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Figure 1. The real part of the bosonic mode w^(y; i i) for t e [0, 10] and K e [0, 4]. 




Figure 3. The real pait of the bosonic mode w^(y; i, i) for t e [0, 20] and K — 0.01. 
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Figure 4. The imaginary part of the bosonic mode W2'(y; 5, 5) for t e [0, 20] and K = 0.01. 
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Figure 5. The real part of the bosonic mode W2~(y; 5, 5) for t e [0, 20] and K = 2. 



Classical Harmonic Oscillator with Dirac-like Parameters and Possible Applications 



13 



0.4 



0.2 - 



-0.2 



-0.4 - 



10 



15 



20 



Figure 6. The real part of the bosonic mode W2'(y; 5, 5) for t e [0, 20] and K 
vertical strip [-0.5, 0.5]. 
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Figure 7. The imaginary part of the bosonic mode Wj" (y; 5 , 5) for t e [0, 20] and K = 2. 
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Figure 9. The fermionic zero mode — 1 / cos t, (red curve), 
(blue curve), for K = 2. 



and the imaginary part of — 1 /w^. 



